f Applications 
in Plant Sciences 



Applications in Plant Sciences 2014 2(2): 1300066 

Primer Note 



Development and characterization of 10 
microsatellite loci in the giant cardon cactus, 
Pachycereus pringlei (Cactaceae)^ 

Carina Gutierrez Flores^, O. Adrian Lozano Garza^, Jose L. Leon de la Luz^, 
AND Francisco J. GarcIa de Leon^'^ 

^Laboratorio de Genetica para la Conservacion, Centro de Investigaciones Biologicas del Noroeste (CIBNOR), Institute 
Politecnico Nacional 195, Col. Playa Palo de Santa Rita Sur, La Paz, Baja California Sur 23096, Mexico 

• Premise of the study: Micro satellite primers were developed for the cardon, a giant columnar cactus {Pachycereus pringlei) 
of the Sonoran Desert, to investigate intraspecific genetic patterns of diversity and population structure. 

• Methods and Results: Using 454 GS-FLX technology and bioinformatics tools, microsatellite primers were successfully identi- 
fied on 282 reads, including di-, tri-, tetra-, penta-, and hexanucleotides. A set of 10 primers were characterized on 80 individu- 
als collected in two areas of the Baja California peninsula, Mexico. All 10 loci were polymorphic, with a mean of 6.3 alleles 
per locus and overall levels of observed and expected heterozygosities ranging from 0.41 to 0.59 and from 0.40 to 0.57, respectively. 
Alleles per individual plant ranged from one to four, suggesting a polyploidal genome. 

• Conclusions: These loci should be useful for future investigations of population structure, genetic diversity, and gene flow in 
the cardon cactus. 
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The giant cardon cactus, Pachycereus pringlei (S. Watson) 
Britton & Rose (Cactaceae), is a columnar cactus that is widely 
distributed in the Sonoran Desert in Mexico, including the Baja 
California peninsula, all islands in the Gulf of California, and 
the coastal region of Sonora State (Turner et al., 1995). This 
species offers an attractive model for evolutionary studies. Its 
wide distribution is characterized by a heterogeneous array of 
environmental conditions (Garcillan et al., 2010), as well as a 
geographic variation in the distribution of its major pollinator, 
the lesser long-nosed bat Leptonycteris curasoae (Fleming 
et al., 1998). Heterogeneity in biotic and abiotic factors can result 
in differential selective pressures across populations, promot- 
ing genetic differentiation and local adaptation. For P. pringlei, 
a latitudinal pattern of abundance, height, and basal diameter in 
the Baja California peninsula is associated with annual and sum- 
mer rainfall, as well as winter temperature (Medel-Narvaez et al., 
2006). This species also exhibits an unusual geographical variation 
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of breeding systems among Cactaceae, including hermaphro- 
ditism (Medel-Narvaez, 2008), gynodioecy (female and her- 
maphrodites), and trioecy (males, females, and hermaphrodites). 
This spatial segregation of sexes is partly associated with varia- 
tions in abundance of its major pollinator (Fleming et al., 1998). 

The Baja California peninsula had a dynamic geologic his- 
tory, including the formation of the Gulf of California approxi- 
mately 5.5 mya (Lonsdale, 1989) and major climatic changes 
during the Pleistocene (Van Devender, 1990) that have had a 
profound effect on the distribution and genetic structure of 
other columnar cacti with similar ranges of distribution, such as 
Lophocereus schottii (Engelm.) Britton & Rose (Nason et al., 
2002) and Stenocereus gummosus (Engelm.) A. C. Gibson & 
K. E. Horak (Clark-Tapia and Molina-Freaner, 2003). Currently, 
there is one molecular study based on low polymorphic markers 
(allozymes) that suggests no genetic differentiation within the 
range of P. pringlei (Fleming et al., 1998). In this study, we 
used the latest 454 GS-FLX titanium sequencing platform to 
locate polymorphic microsatellite loci in P. pringlei and pro- 
pose that these molecular markers be used in genetic studies. 



METHODS AND RESULTS 

Genomic DNA was extracted from fresh cortex tissue of one P. pringlei cac- 
tus growing near La Paz, Baja California Sur, Mexico (23.985°N, 1 10.187°W), 
using a procedure with cetyltrimethylammonium bromide (CTAB) modified 
for mucilaginous cactus tissues (de la Cruz et al., 1997). A voucher specimen 
(JLLL11726 [HCIB28440]) was deposited in the Herbarium HCIB at Centro 
de Investigaciones Biologicas del Noroeste (CIBNOR). Approximately 5 jig of 
genomic DNA was treated with RNAse and submitted to the Genome Sequenc- 
ing Center at the University of California, Santa Cruz, where the sample was 
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Table 1. Characteristics of 10 microsatellite loci developed for the cardon cactus, Pachycereus pringlei. 



Locus 




Primer sequences (5'-30 


Repeat motif 


7;(°c) 


Size range (bp) 


GenBank accession no. 


PpriOl 


F 
R 


GCGAGGGTGTCTCAAATGAT 
TAATCCCCTCCCACCTAACC 


(CT)i2 


58 


202-223 


KC349892 


Ppri02 


F 
R 


TTCCATCGTCCCCTCACTTA 
CATTCACCACCGTGAACACT 


(AG)u 


66 


115-119 


KC349893 


Ppri03 


F 
R 


GGTGTTCCTCGCTCTCATTC 
C T C GC AAAT C C AAGC AAAAT 


(CT)u 


65 


133-155 


KC349894 


PpriU4 


F 
R 


TTCAAATCATTCAATAAAACAAGAGC 
AT GGAGGT C GAGGAT C AAGA 


(TC)i3 


63 


96-112 


Ty^/^O AC\Of\C 


Ppri05 


F 
R 


AAACTGCAGGTGTTTCAGGG 
AATGAAGCGAAAGGAAGCAA 


(GTTT)8 


61 


166-186 


KC349896 


Ppri06 


F 
R 


GCTCACGTTGGCAGATTTGT 
GGTGATGACAAAAGGTTTTGC 


(AAAT)6 


50 


138-146 


KC349897 


Ppri07 


F 
R 


TGGACTTCCAAGGGATAATGA 
TCAACTCAAAGTGTCAGTGCTG 


(AAAT)8 


59 


127-143 


KC349898 


PpriOS 


F 
R 


AATAGCGCATGCCTCAAAGG 

C AAT AGT C C AGAAAT AGGT C AGGT C A 


(CT)8(CTTT)2 


66 


109-126 


KC349899 


Ppri09 


F 
R 


AAGAGAC AGGC C C T GAGAC A 
TCGTAGGTTCCATCACCACA 


(TC)io 


68 


119-140 


KC349900 


PprilO 


F 
R 


TTCTCGAAGCCCCGCTTAC 
GCACGTCAAAGATGCAAGCA 


(GATA)6(GACA)3 


64 


79-107 


KC349901 



Note: = annealing temperature. 



sequenced using 454 pyrosequencing with the GS-FLX titanium platform 
(Roche Diagnostics, Basel, Switzerland), along with seven other species. The 
alignments were run according to the barcode plate strategy, as outlined in the 
GS-FLX Shotgun DNA Library Preparation Method Manual (supplied by 
Roche Diagnostics, Roche Applied Science, Mannheim, Germany). 

All 454 reads were first converted from flowgrams to FASTA and FASTA- 
qual files without 454 adapters, using the sff_extract software (http://bioinf. 
comav.upv.es/sff_extract/index.html). We found 398,286 sequence reads with 
an average read length of 274 bp from two independent runs. MSATCOM- 
MANDER software (Faircloth, 2008) was used to identify perfect microsatel- 
lites containing at least five repeat units for dinucleotides and four for tri-, 
tetra-, penta-, and hexanuclotides. Perfect microsatellite motifs were found in 
10,685 (2.6%) reads, of which dinucleotide repeats were most abundant (62%), 
followed by trinucleotide repeats (31%), and tetranucleotide repeats (6%). 
Using Primer3 software (Rozen and Skaletsky, 2000) embedded in the QDD 
software (Meglecz et al., 2010), primer design was possible in 282 reads, with 
the following parameters: (1) melting temperatures 50-70°C, (2) PGR product 
between 90 and 320 bp, (3) GC content >40%, (4) primer length between 17 
and 27 nucleotides, (5) primer self-complementarities and complement between 
them fulfill the quality criteria used as default parameters in the Primer3 
software. 



Favoring larger motif (repetition numbers) and lower melting temperature 
differences between pairs of primers, we selected 38 loci for PCR amplification 
and a banding pattern of 10 cacti samples collected at El Comitan near La Paz 
(23.985°N, 1 10.187°W). Total DNA was extracted from cortex tissue, using the 
modified CTAB procedure mentioned earlier. PCR reactions were performed in 
10-pL volumes using a thermocycler (ICycler, Bio-Rad Laboratories, Hercules, 
California, USA). Each reaction contained 50 ng of DNA, 0.2 mM of each 
dNTP (no. 18427-088, Invitrogen, Carlsbad, California, USA), 0.4 |iM of each 
primer (Applied Biosystems, Carlsbad, California, USA), Ix PCR buffer 
(20 mM Tris-HCl, 50 mM KCl [pH 8.4]), 1-2.5 mM MgCl2, 0.4 U Tag DNA 
polymerase (no. 1 1615-010, Invitrogen), and 0.1 mg mL"^ of bovine serum albu- 
min (BSA; no. B90015, New England Biolabs, Ipswich, Massachusetts, USA). 
CycHng conditions included initial denaturation at 94°C for 5 min, 30 cycles at 
94°C for 1 min each, 45 s at the locus-specific annealing temperature (see Table 1 
and Appendix SI), extension at 72°C for 1 min, and a final extension at 72°C for 
25 min. PCR products were visualized by silver staining in 6% polyacrylamide 
(Bio-Rad Laboratories), 7.5 M urea, and a 10-bp ladder (Invitrogen). Results of 
PCR amplification of the 38 pairs of primers are given in Appendix SI. 

Ten of the 38 primer pairs yielded distinct bands in the expected size range; 
consistent amplification patterns and fewer stutter bands were selected for char- 
acterization in 40 individuals that were sampled at El Comitan and 40 that were 



Table 2. Genetic diversity of the 10 microsatellite loci developed from the cardon cactus, Pachycereus pringlei.^ 



Locus 




El Comitan {n 


= 40) 






Punta Prieta (n 


= 40) 




Both areas (n = 80) 


A 






He 


A 


Ai 


Ho 


He 


A 


PpriOl 


1 


1-4 


0.67 


0.68 


3 


1-3 


0.38 


0.35 


9 


Ppri02 


2 


1-2 


0.13 


0.39 


3 


1-3 


0.85 


0.63 


3 


Ppri03 


10 


l-A 


0.78 


0.77 


5 


1-4 


0.67 


0.72 


11 


Ppri04 


9 


l-A 


0.77 


0.77 


7 


1-4 


0.74 


0.68 


9 


Ppri05 


5 


1-3 


0.43 


0.57 


3 


1-3 


0.31 


0.35 


6 


Ppri06 


3 


1-2 


0.73 


0.45 


2 


1-2 


0.03 


0.03 


3 


Ppri07 


4 


1-2 


0.75 


0.48 


3 


1-3 


0.68 


0.58 


5 


Ppri08 


4 


1-2 


0.35 


0.37 


2 


1-2 


0.03 


0.03 


5 


Ppri09 


5 


1-3 


0.51 


0.61 


3 


1-2 


0.05 


0.05 


5 


PprilO 


5 


1-3 


0.74 


0.65 


6 


1-4 


0.40 


0.54 


7 


Mean 


5.5 




0.59 


0.57 


3.7 




0.41 


0.40 


6.3 



Note: A = number of alleles per locus; Ai = number of alleles per locus in an individual plant; = expected heterozygosity; = observed heterozygosity; 
n = number of individuals sampled. 

^Locality information for populations: El Comitan, Baja California Sur, Mexico (23.985°N, 110.187°W); Punta Prieta, Ensenada, Baja California, 
Mexico (29.063°N, 114.155°W). 
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sampled at Punta Prieta, Ensenada, Baja California (29.063°N, 114.155°W). 
Duplicate samples of each individual were stored in 96% ethanol at the Labo- 
ratorio de Genetica para la Conservacion at CIBNOR. PGR reactions were 
performed, as described above, using optimized primer-specific annealing 
temperatures (Table 1), magnesium chloride, and forward-labeled primers 
(5' end-labeled with FAM, NED, PET, and VIC; Applied Biosystems). Individu- 
als were genotyped by assessing allele size on a genetic analyzer (model 310, 
Applied Biosystems), using the GeneScan 600 LIZ Size Standard (Applied 
Biosystems) and GeneMapper 4.1 software (Applied Biosystems). 

Locus designation, primer sequences, repeat motifs, annealing temperatures, 
PGR product sizes, and GenBank accession numbers are listed in Table 1. All 
10 loci were polymorphic, exhibiting two to 10 alleles per locus (mean A = 5.5) 
among cacti from El Comitan, and two to seven alleles per locus (mean A = 3.7) 
among cacti from Punta Prieta. Alleles per individual plant ranged from one to 
four in both populations, suggesting a polyploidal genome (Table 2). 

Observed heterozygosity was estimated in Microsoft Excel (Microsoft Cor- 
poration, Redmond, Washington, USA). Expected heterozygosity was esti- 
mated in ATETRA software for tetraploid species that account for all possible 
combinations of allele copy numbers in populations with partial heterozygotes 
(Van Puyvelde et al., 2010). Genetic diversity in terms of heterozygosity varied 
from low to high across loci in both populations. Observed heterozygosity 
ranged from 0.13 to 0.78 (mean = 0.59) and expected heterozygosity ranged 
from 0.37 to 0.77 (mean = 0.57) at El Comitan. At Punta Prieta, observed 
heterozygosity ranged from 0.03 to 0.85 (mean = 0.41) and expected 
heterozygosity ranged from 0.03 to 0.72 (mean = 0.40) (Table 2). For this 
study, we assumed tetraploidy in both populations because we observed as 
many as four alleles per locus in individual plants. Our results are also consis- 
tent with a tetrasomic pattern of inheritance previously suggested in P. pringlei, 
based on an allozymatic study (Murawski et al., 1994) and a chromosome 
counting study (Pinkava et al., 1973). 



CONCLUSIONS 

The 10 microsatellites are the first to be developed for the 
cardon cactus, P. pringlei, and the genus Pachycereus. These 
polymorphic loci will be useful in studies of genetic diversity 
and genetic population differentiation and provide valuable in- 
formation to understand the importance of biotic and abiotic 
factors involved in geographic variation in the genetic pool of 
this widespread cactus of the Sonoran Desert. 
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